function lab4

tspan = [0 10];
load x0.txt;

opcije = odeset('InitialStep', 0.01);

[t1 y1] = ode45(@sustav, tspan, x0, opcije);
plot(t1, y1);
title('Runge-Kutta');

pause;

[t2 y2] = ode23t(@sustav, tspan, x0, opcije);
plot(t2, y2);
title('Trapez');


function s=sustav(t, x)
s = [0*x(1) + 1*x(2)
    -1*x(1) + 0*x(2)];
